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0^ , Abstract. We study the molecular gas content and distribution in the Coma cluster spiral galaxy NGC 4848 which 

\l . is Hi deficient and where the remaining atomic gas is on one side of the galaxy, presumably because of the strong 

>— ^ ' ram-pressure in the rich Coma cluster. Plateau de Bure interferometric CO (1-0) observations reveal a lopsided 

H2 distribution with an off-center secondary maximum coincident with the inner part of the Hi. NGC 4848 is not 
at all deficient in molecular gas as it contains Mh2 ~ 4 10^ M© in the central and inner disk regions. 
As predicted by earlier calculations, ram-pressure has little influence on the dense molecular gas in the inner disk, 
■~^ ■ which appears dynamically normal at our 2 kpc resolution. At the interface between the CO and Hi emission 

i-^ ' regions, about 8 kpc NW of the center, however, strong star formation is present as witnessed by Ha and radio 

continuum emission. From the radio and Ha fluxes we estimate the star formation rate in the northern emission 
Q ■ region to be close to 0.3 M0yr~^, nearly that of an average quiescent spiral. This is the region in which earlier 

5— ( ' Fabry-Perot observations revealed a double-peaked Ha line, indicating gas at two different velocities at the same 

C/3 , sky position. 

C^ ' In order to understand these observations, we present the results of numerical simulations of the galaxy-cluster 

ICM interaction. We suggest that NGC 4848 already passed through the center of the cluster about 4 10^ years 
ago. At the observed stage ram pressure has no more direct dynamical influence on the galaxy's ISM. We observe 
^\. ' the galaxy when a fraction of the stripped gas is falling back onto the galaxy. Ram pressure is thus a short-lived 

JH I event with longer-lasting consequences. The combination of ram-pressure and rotation results in gas at different 

velocities colliding where the double-peaked Ha line is observed. Ram-pressure can thus result, after re-accretion, 
in displaced molecular gas without the H2 itself being pushed efficiently by the intracluster medium. This process, 
however, requires strong stripping, such that only galaxies with radial orbits can be affected as much as NGC 4848. 
A scenario where two interactions take place simultaneously is also consistent with the available data but less 
probable on the basis of our numerical simulations. 

Key words. Galaxies: individual: NGC 4848 - Galaxies: interactions - Galaxies: ISM - Galaxies: kinematics and 
dynamics 

1. Introduction Hi deficient galaxies is comparable to that of correspond- 
ing field galaxies (Kenney & Young 1989 for the Virgo 
Single dish Hi 21 cm line observations show that the spi- ^.j^g^gj.. ^asoh et al. 1991, Boselh et al. 1997 for the Coma 
ral galaxies in the central part (< 1.5°) of a large galaxy cluster). This suggests that the mechanism which is re- 
cluster have less atomic gas than field galaxies of the same sponsible for the removal of the atomic gas does not affect 
morphological type and optical diameter, i.e. they are ^j^g molecular nhase 
Hi deficient (see e.g. Giovanelli & Haynes 1985, Gavazzi 

1989). High resolution observations reveal distortions in The only type of galaxy-cluster interaction which re- 
their Hi distribution (Bravo- Alfaro et al. 2000). Despite moves the neutral Hi gas from the outer disk without mod- 
their reduced atomic gas content, the molecular mass of ifying the molecular gas phase or stellar distribution is 
ram pressure stripping (Gunn & Gott 1972). Ram pres- 
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and is proportional to the galaxy velocity squared and the 
ICM density. Spiral galaxies are recent cluster members 
(TuUy & Shaya 1984, Solanes 2000). Many of them should 
have eccentric orbits, leading them near the cluster cen- 
ter. The gas in clusters is so centrally concentrated that 
ram-pressure is only important close to the center. When 
the consequences of ram-pressure are visible, it means the 
galaxy has already passed through the center if it is Hi 
deficient {DEF > 0.3) (VoUmer et al. 2001). Combes et 
al. (1988) and Kenney & Young (1989) have shown that 
molecular clouds are located too deep in the galaxy's grav- 
itational potential and have column densities too high to 
be removed by ram pressure. Ram pressure stripping re- 
moves only atomic gas clouds with lower column densities 
but it can alter the orbits of the molecular clouds once 
much of the Hi is stripped (VoUmer et al. 2001). 

The effects of ram pressure depend on the galaxy orbit 
and the inclination angle i between the orbital and the disk 
plane. Numerical simulations (VoUmer et al. 2001) show 
that in the case of small i ram pressure not only strips 
the atomic gas from the outer parts of the disk but also 
pushes a considerable amount of gas to smaller galactic 
radii. This leads to a very short-lived enhancement of the 
gas surface density in the inner disk which might lead to 
an enhanced star formation activity there. Furthermore, 
gas clouds which are not accelerated to the escape veloc- 
ity during the stripping event will fall back to the galaxy 
within a few 10® yr, colliding with the galaxy ISM. While 
ram-pressure leads to a long-term decrease in the gas mass 
(Kenney & Young 1989) and star formation rate (SFR; 
Kennicutt 1983), it can temporarily lead to an increase 
in the gas density, and thus presumably the SFR, in two 
ways: (a) for ^ 10^ yr during the closest passage to the 
cluster center or (6) for a few 10* yr after the closest pas- 
sage if re-accretion takes place in disk regions where the 
gas has not been stripped. We have found the signature 
of the back-falling gas in one spiral galaxy in the Virgo 
cluster (NGC 4522, VoUmer et al. 2000). 

The Coma cluster is the densest (tt-co ^ 6 10^^ cm^^, 
Hughes et al. 1989) and most X-ray luminous cluster 
(Lx ~ 3 lO''^^ ergs~^, Jones & Forman 1984) in our neigh- 
borhood. It has been extensively studied at several wave- 
lengths (Hi 21 cm: Giovanelli & Haynes 1985, Gavazzi 
1987;1989, Bravo- Alfaro et al. 2000; CO: CasoU et al. 1991; 
optical: CoUess & Dunn 1996, Biviano et al. 1996; Ha: 
Amram et al. 1992; UV: Donas et al. 1995, X-rays: Briel 
et al. 1992, White et al. 1993, Vikhlinin et al. 1997) re- 
sulting in a detailed view of this cluster. 

In general it is believed that ram pressure strip- 
ping quenches star formation (see e.g. Poggianti 1999). 
However, Donas et al. (1995) observed an enhancement of 
the median UV flux and the fraction of blue star-forming 
galaxies in a ring-like region ^^25' (^0.7 Mpc) from the 
Coma cluster center. The majority of these galaxies are 
Hi deficient. They concluded that a global physical pro- 
cess which leads to the Hi deficiency of these cluster galax- 
ies might induce star formation in a rich cluster such as 
Coma. 



The evolution of the SFR is a key factor determin- 
ing the future appearance of galaxies. Since stars form in 
molecular clouds, we investigate the effect of ram-pressure 
on the molecular gas in the environment where it should be 
most clearly discernible - that of a rich cluster. We present 
high- resolution CO (1-0) observations of NGC 4848, an 
Hi deficient UV-bright spiral galaxy in the richest nearby 
cluster - Coma. It is classified as a blue disk galaxy in the 
sample of Bothun & Dressier (1986). Since it is far away 
from the cluster center (^^1 Mpc) and highly Hi deficient 
but still forms stars at a high rate (^^3 M© yr^^), it is an 
ideal candidate for the scenario described above where the 
stripped gas falls back onto the galactic disc. We are thus 
able to test our simulations directly with observations at 
multiple wavelengths. 

We compare the molecular gas distribution and veloc- 
ity field of NGC 4848 with Hi and Ha interferometric ob- 
servations. With the help of our dynamical model simulat- 
ing ram pressure stripping, we investigate the link between 
the stripping event and star formation. The structure of 
this article is the following: The CO(l-O) observations are 
presented in Sect. 2-4. They are compared in Sect. 5 to 
H band, B band. Hi, Ha, and 20 cm continuum images. 
In Sect. 6 we discuss the observational data. The numeri- 
cal model is presented and compared to the observations 
in Sect. 7. Alternative scenarios are discussed in Sect. 8, 
which is followed by the conclusions (Sect. 9). 

We adopt a distance of 100 Mpc for the Coma 
cluster, corresponding to a Hubble constant of about 
70 kms^^ Mpc^^. All velocities are heliocentric {v = cz). 



2. Observations 

The observations were made with the IRAM interferom- 
eter of Plateau de Bure (see description in Guilloteau et 
al. 1992). We observed NGC 4848 in April 1999 with the 
5D configuration at 112.5674 GHz during 6 hours and in 
May 1999 with the 5C configuration during 7 hours. The 
antenna half-power primary beam was 45"; the primary 
beam field was centered on the optical center of NGC 4848. 
A total bandwidth of 435 MHz (1160 kms~^) centered on 
Vsys=7200 kms~^ was observed with a resolution of 2.5 
MHz (6.6 kms-i). 

Data calibration was made in the standard way us- 
ing the "CLIC" software package (Lucas 1992). The RF 
passband was calibrated at the start of each session ob- 
serving 3C273 or 0923-1-392. The relative phase of the an- 
tennas was checked every 20 mn on the nearby quasars 
1308+326 and 1328+307. The rms atmospheric phase fiuc- 
tuations were typically between 10° and 30°. We applied 
a Banning smoothing to the channels giving 58 channels 
of width 20 kms"'^. The resulting beam size is 5.5" x 3.6". 
The rms noise of one of these channels is 2.3 mJy/beam 
(ATb=0.012 K). 
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Table 1. Physical Parameters of NGC 4848 
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3. Results 




The physical properties of NGC 4848 can be found in 
Table |i|. The complete 12CO(1-0) cube is shown in Fig. |l|. 
The integrated spectrum over a box of 20" x 25" covering 
the whole CO emission can be seen in Fig. |[ The veloci- 
ties are centered on Vsys=7200 kms~^. The flux density of 
the receding part of the galaxy (south) is slightly higher 
than that of the approaching side (north). The total CO 
flux is 5*00=41 Jykms""'^. Using the formula of Kenney & 
Young (1989): 



^ 0.04 - 



A/h, = 1.110* xD^ X 5, 



CO 



(1) 



where D is the distance of the galaxy in Mpc gives a 
total mass of molecular gas of AfH2=4.510^ M©. This 
corresponds to a conversion factor of A^(H2)//co(i-o) = 
2.81020 cm-2K-ikm/s"\ The CO emission distribu- 
tion calculated using a cutoff of 3cr is shown in Fig. ||. 
The major part of the integrated intensity is found in the 
central part of the galaxy. The galaxy's rotation is clearly 



-200 -100 100 200 

heliocentric velocity (km/s) 



Fig. 2. The integrated spectrum over a box of 20" x 25" 
covering the whole CO emission. The velocities are cen- 
tered on vo=7200 kms^^. 



observed in the channel maps. It ranges from -203 to 316 
kms~^ (with respect to vq = 7200 kms~^). The maximum 
of the CO emission coincides with the optical galaxy cen- 
ter. The emission is elongated in the south-east direction 
corresponding roughly to the position angle of the galaxy. 
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Fig. 3. ^^CO(l-O) integrated intensity contours of 
NGC 4848. Contours are (1, 2, 4, 6, 8, 10, 20, 30) x 0.22 
Jy/beamx kms^^. Tfie beam is indicated in the lower left 
corner. 



„ 30 55 



o 

I- 
< 



o 

UJ 

n 




41.0 40.5 40.0 39.5 

RIGHT ASCENSION (B1950) 



Fig. 4. ^^CO(l-O) velocity field. Contours are in steps of 
30 kms~^. The velocities are centered on vo=7200 kms~^. 
The beam is indicated in the lower left corner. 



The most interesting feature is located in the north-west 
of the galaxy center. There, an emission region which is 
hardly spatially resolved is detected. The CO flux within 
this region is 6*00 = 3 Jy km s~^ corresponding to a molec- 
ular hydrogen mass of Mb_2 — 3.3 10^ Mq. This emission 
region appears in nine different channels from -243 to -84 
kms^^. In the integrated intensity map its highest contour 
is at 0.88 Jy/beamxkms~^. 

4. Kinematics 

The velocity field of the CO emission is plotted in Fig. Q. It 
is worth noting that within the northern emission region 
there is a velocity gradient of Av ^ 150 kms~^ within 
only one beam. This gradient appears to be perpendic- 
ular to that of the rotation of the gas in the disk, but 
since this is observed within one beamsize it is not clear 
if the direction of the gradient is real. In order to deter- 
mine the rotation curve we have taken the optical galaxy 
center as a first guess for the dynamical center. First, we 
averaged the approaching and receding side of the galaxy 
up to a radius of 30" excluding points within a sector of 
±30° around the minor axis. The derived dynamical cen- 
ter coincides with the optical one. The results for the other 
parameters are given in Table 0. The parameters as well as 
the rotation curves itself are in good agreement with the 
values derived by Amram et al. (1992). The correspond- 
ing position-velocity diagram together with the our CO 
rotation curve and the Ha rotation curve of Amram et al. 
(1992) are shown in Fig. ||. The rotation curve rises steeply 
in the inner part of the galaxy. In the outer regions it is 
flat with a rotation velocity of tirot=270 kms~^. The ro- 
tation curve rises more slowly at the receding side, which 
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Fig. 5. The p-V diagram along the major axis together 
with the CO and Ha rotation curves. The open boxes 
and circles represent the values obtained by Amram et al. 
(1992). The filled boxes and circles represent our ^^C0(1- 
0) data. The bars indicate the errors on the measurements. 



could be due to the deviation of the galaxy's dynamical 
major axis from the galaxy's position angle and/or due to 
beam smearing. 

5. Comparison with other wavelengths 

5.1. Near Infrared (H band) 

Gavazzi et al. (1996) observed NGC 4848 in the H band 
(1.65 /im) with a resolution of ~ 2" . Fig. g shows an almost 
unperturbed large scale disk, i.e. the isophots are regular 
and do not show a major perturbation. Since the H band 
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Fig. 7. Greyscales: B band CCD image of NGC 4848 with 
a resolution of 0.5"/pixel. Contours: CO emission, for con- 
tour levels see Fig. H. The CO beam is indicated in the 
lower left corner. 



Fig. 6. H band image of NGC 4848 with a resolution of 
~ 2". The contours are from 17 to 22 magarcsec^^ in 
steps of 0.5 magarcsec"^. The arrows indicate four dis- 
tinct spherical emission regions. 



traces the old stellar population, we can conclude that the 
galaxy has not undergone a recent major gravitational 
interaction. However, its gas distribution and kinemat- 
ics traced by the CO, Hi, and Ha emission distributions 
and velocity fields show major asymmetries and pertur- 
bations in the north. This suggests that an ICM-ISM in- 
teraction has selectively disturbed the ISM of NGC 4848. 
Nevertheless, there could be a minor perturbation of the 
old stellar population in the northern part of the galaxy 
where an alignment of four distinct spherical emission re- 
gions, which runs at an angle of ~ 30° across the main 
stellar disk, can be distinguished (arrows in Fig. ^. This 
feature will be discussed in Section W% 



5.2. Optical B band image 

Gavazzi et al. (1990) obtained a CCD image of NGC 4848 
with a resolution of 0.5" in the B band. This image is 
shown in Fig. R together with our CO map. The B band 
image shows a blue galaxy center and a ring-like structure 
which is not centered on the galaxy center but shifted 
~7" to the north. There are four bright blue spots in the 
southern part of this ring and one bright blue spot in 
the north. The southern part of the ring is located at the 
outer edge of the southern elongation observed in CO. The 
bright blue spot in the north lies at the northern edge of 
the detached CO emission region. Moreover, there are two 



faint blue emission regions located at the periphery of the 
CO emission region and another one far away to the north. 

5.3. Hi 

Bravo-Alfaro et al. (1999, 2000) observed several fields 
in the Coma cluster in the Hi 21 cm line with the VLA 
C array. We calibrated and treated the field containing 
NGC 4848 for a second time in order to obtain the maxi- 
mum resolution. The continuum was removed using all ve- 
locity channels between 6800 kms^^ and 7800 kms~^. We 
obtained a final rms noise in one 20 kms^^ wide channel 
of 0.3 mJy per 19" beam. Fig. ^ shows the Hi distribution 
as greyscales and the CO emission as contours. As already 
observed by Bravo-Alfaro et al. (2000) the Hi emission is 
located in the north. With an Hi beamsize of ~ 20" it is 
not spatially resolved. The Hi emission region is located 
at a larger distance from the galaxy center than the CO 
emission region. Since this emission region appears in 5 
different velocity channels above the 2ct level, we are con- 
fident that this detection is real. It represents the neutral 
gas with the highest column density. The lower resolution 
map obtained by Bravo-Alfaro et al. (2000) shows that 
the emission extends northwards. In Fig. we compare 
the spectra of the Hi and CO emission regions integrated 
over the area where emission is detected. Both spectra 
show emission in the same velocity range between 7000 
and 7100 kms^"'^, the CO line being two times larger. This 
confirms that both detections are real. Furthermore, the 
integrated Hi spectrum of Giovanelli & Haynes (1985) has 
the same shape and width as our CO spectrum. Thus, the 
northern CO emission region and the Hi emission region 
are very probably associated. Given the small errors on 
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Fig. 8. Greyscales: Hi distribution of 19" resolution. Fig. 10. Greyscales: Ha distribution with a resolution of 
Contours: CO emission, for contour levels see Fig. ^. The 2". Contours: CO emission, for contour levels see Fig. ^ 
CO beam is indicated in the lower left corner. The CO beam is indicated in the lower left corner. 
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Fig. 9. Solid line: normalized CO integrated spectrum of 
the northern emission region. Dashed line: normalized Hi 
integrated spectrum of the northern emission region. 



the optical and Hi velocities (^20 kms^^) compared to 
their difference (Aw ~ 180 kms~^, see Tab. |^), we think 
that the offset between the optical and Hi velocities is real. 
The Hi flux in the C-array VLA map of this region 
is S'hi=0.05 Jykms-i (Mhi=1.210* Mq), whereas the 
single-dish Arecibo observations show that much more low 
surface density Hi is present (Mhi=1.7109 Mq). The CO 
flux of the northern emission region is 5*00=3 Jykms"^ 
(Mh2=3.3 10* M©), which represents -7% of the total CO 
flux. More sensitive Hi observations are necessary to de- 
termine whether the gas is chiefly atomic or molecular in 
this region. 



5.4. Ha 

We compare here our CO data with the Ha Fabry-Perot 
observations of Amram et al. (1992). Their map is shown 
as greyscales together with the CO data as contours in 
Fig. |l^. We determined the coordinates of the galaxy cen- 
ter of the Ha image by fitting it to the B band image. The 
Ha emission distribution shows a ring structure very sim- 
ilar to the one observed in the B band image. However, it 
contains only two bright spots in the south west and in the 
north where a lot of massive stars are formed. The other 
three bright spots in the south of the B band image do not 
have counterparts in Ha. The central hole is probably due 
to a too severe continuum subtraction in this region as a 
comparison with the Ha image of Gavazzi et al. (1998) in- 
dicates. The northern bright spot is located near the CO 
emission region. The uncertainty in the determination of 
the galaxy center coordinates in the Ha image does not 
allow us to conclude if both emission regions coincide ex- 
actly. We conclude here that both are located very close 
to each other within the same region. The northern Ha 
emission region contains ~20% of the total Ha emission. 
In the south, the outer edge of the CO elongation nicely 
fits the Ha distribution. 



5.5. Ha /Ha continuum 

Gavazzi et al. (1998) obtained a narrow band Ha CCD 
image of NGC 4848 with a resolution of 1.2". We show 
in Fig. |l^ the ratio between the Ha line emission and its 
continuum, in order to study the nature of the off-plane 
NIR emission regions (Fig. H). The image shows clearly 
the ring of highest Ha emission (see Fig. |0[). Moreover, 
the western part of the galaxy next to the Ha ring has 
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Fig. 12. Greyscales: 20 cm continuum distribution of 5" 
resolution. Contours: CO emission, for contour levels see 
Fig. H. The CO beam is indicated in the lower left corner. 

figuration (FIRST; Becker et al. 1994). These maps have 
1.8" pixels, a typical rms of 0.15 mJy, and a resolution of 
5". The 20 cm continuum map of NGC 4848 together with 
the CO emission map is shown in Fig. |l2|. The 20 cm con- 
tinuum map shows enhanced emission to the north west. 
It coincides exactly with the maximum of the CO emission 
in the northern emission region. The southern elongation 
of the CO emission does not have a counterpart in the 20 
cm continuum map. 



Fig. 11. Ratio between the Ha line emission and its con- 
tinuum. The contours are from 0.2 to 4 in steps of 0.2. We 
have cut the south-eastern part of the galaxy because of 
straylight problems due to a nearby star (Gavazzi et al. 
1998). 

an enhanced Ha/continuum ratio. The NIR emission re- 
gions, which are running at an angle of ~ 30° across the 
galaxy's disk (see Fig. ^), have a low Ha/continuum ra- 
tio, similar to that of the extended disk. The north west- 
ern Ha emission maximum is not aligned with the NIR 
emission regions but 2.5" ~ 1 kpc offset to the south. On 
the contrary, in the scenario of a tidally stretched dwarf 
whose ISM is interacting with the ISM of NGC 4848, one 
would expect (i) an alignment of all these regions and (ii) 
a higher Ha/continuum ratio in one of the two brightest 
NIR off-plane regions (see, e.g., Clemens et al. 2000). 

5.6. 20 cm 

The Very Large Array (VLA) carries out a systematic sur- 
vey of the northern sky at 20 cm wavelength in the B con- 



6. Discussion of the observations 

Many of the physical properties of NGC 4848 are compa- 
rable to those of Virgo cluster Sab galaxies with the same 
blue luminosity and the same Hi deficiency (Kenney & 
Young 1989): 

- S(60 /zm)/S(100Aim)=0.52 

— ratio of molecular hydrogen mass to optical area: 
log(MH,/i?2^t)=0.37Mopc-2 

— ratio of total gas mass to optical area: log((AfHi + 
MHj/^^pt)=0.51Mopc-2 

- ratio of Hi to H2 fiux: log(S'co/S'Hi)=l-8 , 

whereas its Ha line fiux is a factor of at least two higher 
than that of a Virgo cluster spiral of the same morpholog- 
ical type, the same diameter, and the same Hi deficiency. 
In addition, its equivalent width is about 3 times higher 
than that of "normal" field Sab galaxies (Kennicutt 1983) 
and resembles that of an Sc spiral with a star formation 
rate of 



SFR 



i(H„) 



1.26104iergs- 



4.5 10-'''^ L(FIR) ~ 3 Mq yr-i(2) 



(Kennicutt et al. 1994, Kennicutt 1998). 



VoUmer et al.; ^^ CO (1-0) observations of NGC 4848 



Calculating the CO deficiency 



Def 



CO 



l0gAf(i/2)c-l0gAf(i/2)o 



(3) 



where M{H2)o is the observed and M{H2)c the expected 
molecular mass (see Boselli et al. 1997) gives negative val- 
ues, i.e. it is not CO deficient. 

NGC 4848 is a Hi deficient galaxy (Defni ^ 0.5; 
Gavazzi 1989; Giovanelli & Haynes 1985), i.e. it has lost 
about 2/3 of its atomic gas. The fact that the disk to 
bulge ratio in the H band appears to be small (Fig. 0) 
also indicates that NGC 4848 was likely an Sc galaxy be- 
fore entering the cluster. Since there is no identified major 
optical companion neither a major asymmetry in the H 
band image, a major gravitational interaction can be ex- 
cluded. The alignment of spherical emission regions in the 
NIR and all visible wavelengths might represent a tidally 
sheared dwarf galaxy which crosses the disk of NGC 4848. 
If this dwarf galaxy has been gas rich this could have an 
influence on the gas distribution of NGC 4848. But since 
the dwarf galaxy is very small compared to NGC 4848 
(Fig. H), it is highly improbable that it is responsible for 
the displacement of more than 10^ Mq of atomic gas to 
the north of NGC 4848. Therefore, ram pressure stripping 
is the most probable cause for the atomic gas removal. 
We will discuss the scenario of a tidally stretched dwarf 
further in Sect. B.2. 



Numerical simulations (VoUmer et al. 2001) indicate 
that Hi deficient galaxies (Dcfni > 0.3) which are strongly 
affected by ram pressure stripping and which are located 
far away from the cluster center are observable most 
probably after their closest passage to the cluster cen- 
ter. Ram pressure Pram = PiCMfg^i is only high enough 
near the cluster center, where the ICM density picM 
and the galaxy's velocity Ugai increase rapidly, to over- 
come the restoring gravitational forces. Thus, ram pres- 
sure only acts in the immediate vicinity of the cluster 
center. In addition, it needs only several lO'' yr to ac- 
celerate the Hi clouds to their escape velocity (Murakami 
& Babul 1999), thus stripping happens very rapidly. Since 
a galaxy passes the cluster center with a very high ve- 
locity (Ugai > 2000 kms~^), the probability of observing 
it when stripping sets in is very low. It is therefore most 
likely that we observe an Hi deficient galaxy, which has a 
distorted Hi distribution, when it comes out of the cluster 
center. We will adopt this view in this work. Furthermore, 
Vollmer et al. (2001) showed that there is re-accretion of 
stripped material, which has not been accelerated to the 
escape velocity, between 2 and 5 10® yr after the closest 
passage. This infalling gas hits the galactic disk leading 
to frequent cloud-cloud collisions. The compression of the 
ISM due to these collisions might lead to an enhancement 
of the local star formation activity where the accreted gas 
falls on the disk. 

NGC 4848 is located at 26'(~0.75 Mpc) from the clus- 
ter center. Assuming a mean velocity of 1800 kms~^, the 
closest passage was about 4 10® yr ago. Therefore we as- 
sume that the galaxy is in the stage where material is 
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Fig. 13. Spectra of the Ha emission region in the north of 
the galaxy; x-axis: velocity (6940 to 7468 kms^^); y-axis: 
Ha emission intensity in arbitrary units (data cube from 
Amram et al. 1992). 



falling back on its disk. The shocks produced by the en- 
counters of infalling clouds with the ISM will compress the 
gas and star formation will take place. Given this picture, 
we can interpret the observational data shown above. 

The Hi and CO maps show emission in the north of 
the galaxy. This means that the galaxy's gas content of 
highest column density in the disk is located in the north 
of the galaxy center. We propose that the compression 
leads to a phase transition from atomic to molecular gas 
in the region where the infalling gas hits the disk, i.e. in 
the north. 

The shift of the Ha emission ring to the north indi- 
cates an enhanced massive star formation in the north of 
the galaxy. However, the extended old disk component is 
symmetric (Gavazzi et al. 1990). The bulk of the neutral 
gas moves around the galaxy center with an almost con- 
stant rotation curve while the ionized gas content shows 
an important deviation: in the north western part of the 
disk the Ha emission has a double peaked profile (Fig. 03 
extracted from the Ha data cube of Amram et al. 1992). 
The peak at lower velocities (left) is due to the regular ro- 
tation of the gas within the galaxy following the rotation 
curve (Fig. ||). The stronger peak (right) at higher ve- 
locities is separated by ~-|-100 kms~^ from the first one, 
which corresponds to the disk rotation. The Ha line flux 
in this region represents ~ 20% of the total Ha line flux. 

We interpret this feature as the consequence of in- 
falling gas which collides with the ISM within the galaxy. 
The cloud-cloud collisions lead to the compression of the 
gas and the magnetic field. This gives rise to an enhanced 
star formation observed in the Ha map. The enhanced 
star formation activity together with the compressed mag- 
netic field leads to the northern enhancement in the 20 cm 
continuum map. The enhanced number of young and hot 
0/B stars enhanced the ambient UV radiation field in the 
norther part of the galaxy, leading to an increase of the gas 
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temperature. This effect, together with the increase of the 
molecular gas mass due to shock compression during the 
re-accretion, increases the ^^ CO (1-0) emissivity. The line 
emission thus becomes observable in the northern part of 
the outer disk. 

In the next Section we present a detailed dynamical 
simulation of the infalling galaxy including ram pressure 
stripping in order to test the hypothesis made above. 



7. Dynamics 

7.1. The numerical model 

Very few simulations have been done to quantify ram 
pressure stripping using Eulerian hydrodynamic (Takeda, 
Nulsen, & Fabian 1984, Gaetz, Salpeter & Shaviv 1987; 
Balsara, Livio, & O'Dea 1994) or SPH codes (Tosa 1994; 
Abadi, Moore, & Bower 1999). All methods have their 
advantages and their limits. Since we need to model an 
extreme density gradient between the ICM and the ISM 
in three dimensions, one would need a very high number 
of SPH particles or a very high resolution for Eulerian 
grid codes which makes calculations very memory, and 
time, consuming. The sticky particles code we use has 
the advantage that the effects of the ICM are directly 
implemented as an additional acceleration of the clouds 
which are exposed to ram pressure. The gas viscosity is 
due to direct (collisions) or indirect (gravitational) inter- 
actions between the clouds. We are the first to apply a 
time-dependent ram pressure profile, i.e. a Gaussian pro- 
file. The consequences of the decrease of ram pressure after 
the galaxy's closest passage to the cluster center is that 
material which has not been accelerated to the escape ve- 
locity falls back to the galaxy. This has not been observed 
in simulations before, because of the use of a constant, 
time-independent ram pressure. 

We used the three-dimensional N-body code described 
in detail in VoUmer et al. (2001). The particles repre- 
sent gas cloud complexes which are evolving in an ana- 
lytically given gravitational potential of the galaxy. This 
potential consists of two spherical parts: the dark mat- 
ter halo and the stellar bulge, and a disk potential (Allen 
& Santillan 1990). The model parameters are: (i) halo: 
03=12 kpc, A/3 = 2.110" M0, (ii) bulge: 6i=400 pc, 
Ml = 2.8101" Mq, (iii) disk: a2=5 kpc (disk scale 
length), 62=250 pc (height from the disk plane), M2 = 
1.710^1 Mq. The outcoming velocity field has a constant 
rotation curve of Wiot ^^250 kms"^ (Fig- |lj)- The optical 
diameter of NGC 4848 D ^^ 44 kpc corresponds to --9 disk 
scale lengths. 10 000 particles of different masses are ro- 
tating within this gravitational potential. The cloud mass 
distribution is n{m) ex: m^^'^ (Knude 1981) within a range 
of ~ 10* to 510^ Mq. 

A radius which depends on the mass is attributed to 
each particle. During the disk evolution they can have 
inelastic collisions. The outcome of these collisions is sim- 
phfied following Wiegel (1994): 



model rotation curve 



radius (pc) 

Fig. 14. The model rotation curve. The error bars indicate 
the velocity dispersion of the particles. 



— for ri — r2 < fe < ri + 7*2 : 
fragmentation 

— for h < ri — r2 and Wosc > Vf'- 
mass exchange 

— for h < ri — r2 and Vcsc < ''^f : 
coalescence, 

where h is the impact parameter, ri and r2 the cloud radii, 
Wf the final velocity difference, and Wosc is the escape ve- 
locity. This results in an effective disk viscosity. The next 
neighbor search and selfgravitation between the clouds is 
done with a three-dimensional treecode (Barnes & Hut 
1986). 

As the galaxy moves through the ICM its clouds are 
accelerated by ram pressure. Within the galaxy's iner- 
tial system its clouds are exposed to a wind coming from 
the opposite direction of the galaxy's motion through the 
ICM. The effect of ram pressure on the clouds is simu- 
lated by an additional force on the clouds in the wind 
direction. Only clouds which are not protected by other 
clouds against the wind are affected. 

When the galaxy approaches the cluster center, its ve- 
locity increases. At the same time the surrounding ICM 
density increases. This leads to an increase of the ram pres- 
sure on the ISM clouds Pram = PiCM^gaij where picM is the 
ICM density and v^a\ is the velocity of the galaxy. We take 
this evolution of Pram into account in adopting the follow- 
ing profile Pram(i) = 28poWo^xp( - (i/1.210'*)^), where 
po = lO^"' cm~'^ and wo=1000 kms"^, where t = yr is 
the time of the galaxy's closest approach to the cluster 
center. 

With the scenario in mind we have taken a snapshot 
of the evolved galaxy 400 Myr after the closest passage 
to the cluster center. VoUmer et al. (2001) have shown 
that the resulting Hi deficiency increases with increasing 
inclination angle i between the galaxy's disk and its or- 
bital plane, whereas the fraction of re-accreted material 
increases with decreasing i. In our case, a considerable 
fraction of stripped gas clouds are accelerated by ram pres- 
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sure to velocities below the escape velocity. In addition, for 
small i the re-accreted material hits the disk mainly at one 
single point. For increasing i the area where re-accretion 
takes place increases and the gas falls from the z direc- 
tion onto the disk. For the comparison with observations 
we have therefore chosen i=20° and a ram pressure max- 
imum of 28po''^o- This satisfies all necessary conditions: a 
high stripping efficiency, a relatively high re-accreted mass 
fraction and a small area where the re-accreted cloud com- 
plexes hit the disk ISM. The final Hi deficiency of this 
simulation is DefHi=0.7. 

Fig. |l^ shows the evolution of the model galaxy in 
timesteps of Ai '^ 70 Myr. The ram pressure maximum 
occurs at i=0 yr. 

Were NGC 4848 falling into the cluster core, it would 
be even less disturbed than in the first panel of Fig. |l^, 
which is for 98 Myr before reaching the core. This is 
consistent with the conclusion of VoUmer et al. (2001) 
that highly distorted, Hi deficient galaxies which are most 
likely located far away from the cluster center (and in 
particular NGC 4848) are coming out of the Coma cluster 
core. 

At t '^ 100 Myr before the galaxy's closest passage to 
the cluster center a density enhancement grows in the di- 
rection of the galaxy's motion. The gas of the outer Hi disk 
is driven away from the galaxy. Since it still has the angu- 
lar momentum of rotation, it forms a ring-like structure. 
At the closest passage to the cluster center an accelerated 
arm in the north-east and a decelerated arm in the south 
are formed. Both arms rotate counter-clockwise. The ro- 
tation of the southern arm is stopped by the fading ram 
pressure. At t '^ 200 Myr after the galaxy's closest pas- 
sage to the cluster center the stripped gas which has been 
accelerated to velocities remaining below the escape ve- 
locity begins to fall back onto the galaxy following both 
arms. The infalling gas ffows cross in the north-west of the 
galaxy where they form an expanding gas shell. 

The model snapshot at i = 400 Myr after the closest 
passage of the cluster center is shown in Fig. 16, The direc- 
tion of the galaxy's motion indicated in Fig. 15. The hole 
in the center of the galaxy is artificial, because the clouds 
in the very center are arbitrarily removed from the model. 
Since the timestep is roughly inversely proportional to the 



acceleration a, Ai 



i?w^Qt, where R is the galac- 



tic distance and we assume Wrot^const, this procedure in- 
sures a large enough timestep (Ai ~ 10^ yr). The galaxy's 
main Hi disk is truncated at a diameter of Dhi ~ 15 kpc. 
Assuming an optical diameter of 40 kpc for the model 
galaxy, this gives a ratio between the Hi and the optical 
diameter of DniZ-Dopt ~ 0.4, which fits in the correspond- 
ing plot of Cayatte et al. (1994). If we assume the same 
stripping mechanism acting in the Coma cluster as for the 
Virgo cluster, we expect a ratio between the Hi and the 



optical diameter to be Dju/D, 



opt 



0.4 - 0.5. With the 
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optical diameter of NGC 4848 this leads to a predicted 
Hi diameter of /jP^j^dictcd _ 20 kpc (40"), which is only 
slightly larger than the model Hi diameter. The detected 



Fig. 15. Snapshots of the galaxy evolution. The elapsed 
time is indicated at the top of each panel. The position 
and inclination angle of NGC 4848 are used. The galaxy 
rotates clockwise. It is moving in the north-west direction, 
i.e. the wind is coming from the north-west (indicated by 
the arrows). The length of the arrow is proportional to 
ram pressure {pv^). 



Hi emission region is located at a distance from the galaxy 
center of ^ 15", thus well within this predicted diameter. 
A shell of high surface density Hi gas is located in the 
north (Fig. Hq, left panel). It spans an angle of ^ 45°. This 
shell is located mainly above the disk plane (Fig. O, right 
panel). It consists of gas which has been removed from 
the disk by ram pressure in the wind direction (20° off 
plane). This material is now falling back into the galaxy. 
In addition, a tenuous gas component can be found north 
and south of the galaxy. 
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Fig. 16. Model gas distribution. Left: Snapshot of the simulation at i = 400 Myr after the closest passage of the 
cluster center. The sizes of the circles correspond to the sizes of the cloud complexes. The clouds are projected on the 
x-y-plane. Right: The clouds are projected on the x-y-plane. 
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Fig. 17. Model velocity field. For clarity only 10% of the vectors are shown. Left: model velocity field projected on 
the x-y-plane. The galaxy rotates counter-clockwise. The length of the arrows is proportional to the absolute value of 
the projected velocity. Right: model velocity field projected on the a;-z-plane. 



The projection of the cloud velocities on the x-y-plane 
and a;-2:-plane are shown in Fig. 17, The galaxy rotates 
counter-clockwise. The tenuous gas component at large 
distances from the galaxy falls back to the galaxy. The 
northern shell consists of two counter-rotating streamers. 
The major streamer comes from the west in the sense of 



the galactic rotation, the minor one arrives from the east 
in a counter-rotating sense. They correspond to the two re- 
accretion arms which form after the stripping event when 
ram pressure has already ceased. The first arm rotates in 
the sense of the galaxy rotation and meets the disk in 
the south, the second arm is counter-rotating and meets 
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Fig. 18. Comparison between the simulation snapshot 

and the Hi observations of Bravo- Alfaro et al. (2000). 

Greyscale: model column density distribution. Contours: Fig. 19. Ha emission distribution. Greyscale: Model. 

Hi emission distribution of Fig. |. Contours: observations (Fig. |^. 



the galactic disk in the north. The northern shell is thus 
expanding to the north. 

7.2. Comparison with observations 

We have generated the observer's view of this gas distribu- 
tion with PA=158° and j=75°. For the comparison of the 
model column density distribution with the Hi data we 
assume that the gas within a distance of 7 kpc (15") from 
the galaxy center is mainly molecular whereas the gas at 
larger distances is mainly atomic. We can thus compare 
directly both distributions (Fig. |l^). The low resolution 
Hi distribution of Bravo- Alfaro et al. (2000) suggests that 
the bulk of the atomic gas is located more to the north 
than the emission region shown in Fig. g. The model Hi 
emission region is located in the right direction and at 
approximately the right distance from the galaxy center. 
Thus, we suggest that the observed Hi emission regions 
corresponds to the inner part of the expanding northern 
shell of our simulation. 

In order to compare the Ha image with the model, we 
again assume that star formation takes place where the re- 
accreting material hits the disk, i.e. where we have a dou- 
ble line profile. Each cloud complex in the model snapshot 
was checked for a velocity which leads to a collision with a 
neighboring cloud complex. We thus project the velocity 
vectors (vi, V2) of each pair of clouds on their relative 
distance vector y^^'. v^"'^ = vi ■ ri2 ; ^2'°"' ~ '^2 • ri2 ■ 
We apply a weighting factor to each cloud which is equal 



to the number of possible collisions fulfilling the criterion: 
ri2 < 1 kpc, -yf™^ has the opposite sign of ^2™''- All cloud 
complexes which will not collide are removed from the 
cube. At distances R < 7.5 kpc, i.e. within the ring ob- 
served in Ha, the mean free path of a cloud is smaller than 
outside the disk. The criterion ri2 < 1 kpc might thus not 
be valid at small distances. Moreover, in the inner disk, 
star formation is more likely due to density waves within 
the stellar and gaseous disk. Since this aspect is not in- 
cluded in the model we prefer to give the clouds a,t R < 7.5 
kpc a uniform weight (Vollmer et al. 2000). The model 
emission inside the galaxy disk is thus density weighted, 
the emission outside the disk is collision weighted. 

The result of this procedure together with the Ha emis- 
sion distribution can be seen in Fig. |l^. The extent of the 
model and observed Ha emission fits remarkably well. The 
two brightest spots in the model are located in the north- 
west and in the south of the galaxy center. Especially the 
north-western maximum coincides exactly with the ob- 
served maximum of the Ha emission. 

In order to check if we can reproduce the double line 
profile in the north-eastern emission region as it has been 
observed by Amram et al. (1992), we made spectra of the 
model cube in this region (Fig. EO) , which can be directly 
compared with the Ha spectra of Fig. |l^. As in the ob- 
served spectra, we find a double line profile in the model 
cube. The main difference between the model and obser- 
vations lies in the different intensities of the lines. The 
observed line of highest intensity is not due to regular 
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Fig. 20. Model Ha spectra of the region in the west of the 
galaxy, x-axis: velocity (6940 to 7468 kms— 1 same as in 
Fig. 03); y-axis: gas surface density in arbitrary units. 



rotation, whereas the model shows the contrary. This dif- 
ference is not significant, however, because we have not 
included an explicit recipe for star formation or radiative 
transfer in our model, making it impossible to predict the 
outgoing ionizing UV field in the star forming region. 

7.3. Towards a coherent scenario 

Our model includes only 10 000 particles but the calcula- 
tion is very long for each run due to the small time steps, 
limiting the exploration of parameter space. Our aim is 
to show generic features which arise naturally from our 
simulations. These effects are more or less pronounced for 
different i and Pmax- The val ues chosen represent a corn- 
Nevertheless, there are 



promise as described in Sect. 7.1 



typical features which appear only in a given range of i 
and Pmax, as the amount of re- accretion and the resulting 
Hi deficiency. 

We are able to reproduce 

— the observed Hi distribution, i.e. that most of the Hi 
is found in the north of the galaxy, 

— the north-western maximum of the observed Ha emis- 
sion distribution, 

— the observed double Ha line profile in the north-west 
of the galaxy, 

— the enhanced 20 cm radio continuum emission region 
to the north of the galaxy center. 

In our model the re-accreting gas hits the disk at the 
position of the observed CO emission region. We expect 
that the molecular gas has condensed from the colliding 
atomic gas, inducing the observed star formation. The ob- 
served radio continuum peak near the CO/Hi emission re- 
gion arises naturally from an ISM-ISM collision in which 
the magnetic fields are compressed at the same time as 
the clouds. 



8. Alternative scenarios 

The scenario presented here, with the observational asym- 
metries of the gas distributions and perturbations of the 
velocity fields, is not unique. The distance to the Coma 
cluster is such that the resolution of the presented obser- 
vations is lower than that of the model. Therefore, we will 
discuss two alternative scenarios. 



8.1. Ram pressure stripping in action 

If ram pressure is stripping the gas now, we could imagine 
a scenario where NGC 4848 approaches the cluster center. 
The neutral hydrogen and maybe a part of the molecu- 
lar gas (Kritsuk 1983) are displaced to the north west of 
the galaxy by the ram pressure of the ICM. In this way 
NGC 4848 would resemble NGC 4419 in the Virgo clus- 
ter (Kenney et al. 1990), which also shows an asymmetric 
distribution of the molecular gas. Nevertheless, there are 
important differences. NGC 4419 is red {B - V)t = 1.0 
(Kenney et al. 1990), whereas NGC 4848 is a blue disk 
galaxy {B - V)t = 0.44 (RC3). Moreover, the asym- 
metry of the CO distribution is much more pronounced 
for NGC 4419 (a factor of 2 in CO Imninosity) than for 
NGC 4848 (a factor of 0.1). Thus, they might not be di- 
rectly comparable. 

For the scenario of active stripping, two main questions 
remain: (i) The galaxy has lost about 2/3 of its atomic gas. 
Thus ram pressure stripping has already been very effi- 
cient at larger distances from the cluster center (> 1 Mpc) 
in order to have removed more than 3 10^ Mq. This gas 
must have been ionized, because we do not detect it in the 
21 cm line. At these distances the ICM density is very low 
{ric ^ 10^^ cm~^). If we assume an Hi column density of 
Shi = lO^*' cm^^, a rotation velocity Uiot = 270 kms^^, 
and a galactic radius of i?=10 kpc, we obtain from the en- 
ergy density equation for face-on stripping (Gunn & Gott 
1972) 



SgasW.^ot^ ^ = ™ 



p'T-c^gal J 



(4) 



where rrip is the proton mass, and we have a galaxy veloc- 
ity of Ugai ~ 1600 kms^^. The Keplerian velocity for cir- 
cular orbits at 1 Mpc from the cluster center, assuming an 
enclosed total mass of Aftot = 5 lO^'' Mq (Hughes 1989), 
is WKcpicr ^ 1500 kms^^. If the galaxy was stripped face- 
on, the scenario of active ram pressure stripping would be 
consistent with the Hi 21 cm line observations. 

NGC 4848 has a low radial velocity with respect to 
the cluster mean (Au ^ 350 kms~^). Thus, the galaxy's 
velocity vector must have a large component in the direc- 
tion of the cluster center. Since it is highly inclined, it is 
most probable that NGC 4848 is stripped nearly edge- 
on. In this case, Eq. does not apply, i.e. ram pressure 
must be much higher to strip the outer atomic hydrogen. 
According to our model simulations (Fig. 15), a displace- 
ment of the gas is only clearly detectable after the galaxy's 
passage through the cluster center. Even if we underesti- 
mate the ram pressure efficiency, it is unlikely that a dis- 
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placement already takes place at At ^ 400 Myr before the 
galaxy's closest approach to the cluster center. 

8.2. Tidal interaction with a dwarf galaxy 

Fig. y shows an alignment of four nearly spherical emission 
regions at an angle of ~ 30° with respect to the galaxy's 
position angle. This could be a tidally stretched dwarf 
galaxy crossing the disk, in the same way the Sagittarius 
dwarf crosses the Galaxy (see Edelsohn & Elmegreen 1997 
for numerical simulations). In this case, one would expect 
a major emission region which corresponds to the core of 
the dwarf galaxy. The two strongest emission regions at 
the extremities of the alignment have approximately the 
same brightness in the H band, weakening the interpreta- 
tion of the alignment as a tidally stretched dwarf galaxy. 

Even if the galaxy was gas rich, this can not explain the 
displacement of ^ 10^ M0 of atomic hydrogen. Moreover, 
in an ISM-ISM interaction one would expect strong Hi 
tails following the stellar tails (see, e.g., Weliachew et al. 
1978). An ISM-ISM collision could explain the detached 
northern CO emission region together with the enhanced 
radio emission, but it cannot account for the higher Ha 
line to continuum ratio in the west of the Ha ring, either 
for the large spatial extent of the region where Ha double 
line profiles were observed, or for the offset between the 
alignment and the western Ha maximum. 

As ram-pressure stripping is necessary in the two pos- 
sible scenarios and naturally accounts for the observations 
according to our simulations, it is currently the preferable 
explanation. Further observations (e.g., with the HST) 
may allow us to determine whether a small galaxy is also 
being tidally torn apart. 

9. Conclusions 

We have shown ^^CO(l-O) Plateau de Bure observations 
of the highly Hi deficient Coma galaxy NGC 4848. A de- 
tached CO emission region in the north of the galaxy cen- 
ter was detected. This region is located near an Ha emis- 
sion maximum where an Ha double line profile was ob- 
served earlier. Furthermore, an enhancement in the 20 cm 
radio continuum coincides with the detached northern CO 
emission region. 

We are able to reproduce naturally the main charac- 
teristics of the emission distributions and velocity fields 
of the multiple wavelength observations with the help of 
a numerical simulation including ram pressure stripping. 
In this scenario, the galaxy passed near the cluster center 
~400 Myr ago. NGC 4848 is emerging from the cluster 
center with a radial velocity of ^^-400 kms~^ and is lo- 
cated in front of the cluster center. We propose that it was 
an Sc galaxy which entered the cluster for the first time 
~1 Gyr ago. During its recent passage through the cluster 
core most of its atomic gas was removed by ram pressure 
stripping. Its molecular gas has not been stripped because 
it is located deep in the gravitational potential well of the 
galaxy and because of its very high column density. The 



atomic gas which was not accelerated to a velocity above 
the escape velocity now falls back onto the galaxy. The 
compression of the neutral gas due to accretion leads to 
a phase transition from atomic to molecular gas and to 
star formation activity. We expect the star formation rate 
to decrease within the next Gyr due to the lack of a gas 
reservoir in the outer disk. 

An alternative scenario, where ram pressure stripping 
is active and a tidally stretched dwarf galaxy crosses the 
disk at the same time, is also consistent with the avail- 
able data but less probable on the basis of our numerical 
simulations. 
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